Search for: All records

Abstract A spinning black hole (BH) accreting from a disk of strongly magnetized plasma via a magnetically arrested disk is known to produce an efficient electromagnetic jet powered by the BH’s spin energy. We present general relativistic radiative magnetohydrodynamic simulations of magnetically arrested systems covering a range of sub- to super-Eddington accretion rates. Using the numerical results from these simulations, we develop formulae to describe the magnetization, jet efficiency, and spin evolution of an accreting BH as a function of its spin and accretion rate. A BH with near-Eddington accretion experiences a mild degree of spin-down because of angular momentum loss through the jet, leading to an equilibrium spin of 0.8 rather than 1.0 at the Eddington limit. As the accretion rate increases above Eddington, the spin-down effect becomes progressively stronger, ultimately converging on previous predictions based on nonradiative simulations. In particular, spin evolution drives highly super-Eddington systems toward a BH spin near zero. The formulae developed in this letter may be applied to galaxy- and cosmological-scale simulations that include BHs. If magnetically arrested disk accretion is common among supermassive BHs, the present results have broad implications for active galactic nucleus feedback and cosmological spin evolution. 

ABSTRACT We present two general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of magnetically arrested discs (MADs) around non-spinning (a* = 0) and spinning (a* = 0.9) supermassive black holes (BHs). In each simulation, the mass accretion rate is decreased with time such that we sample Eddington-scaled rates over the range $$3 \gtrsim \dot{M}/\dot{M}_{\rm {Edd}}\gtrsim 0.3$$. For the non-spinning BH model, the total and radiative efficiencies increase as the accretion rate decreases, varying over the range $$\eta _{\rm {tot}}\sim 9\!-\!16{{\ \rm per\ cent}}$$ and $$\eta _{\rm {rad}}\sim 6{-}12{{\ \rm per\ cent}}$$, respectively. This model shows very little jet activity. In contrast, the spinning BH model has a strong relativistic jet powered by spin energy extracted from the BH. The jet power declines with accretion rate such that $$\eta _{\rm {jet}}\sim 18{-}39{{\ \rm per\ cent}}$$ while the total and radiative efficiencies are $$\eta _{\rm {tot}}\sim 64{-}100{{\ \rm per\ cent}}$$ and $$\eta _{\rm {rad}}\sim 45{-}79{{\ \rm per\ cent}}$$, respectively. We confirm that mildly sub-Eddington discs can extract substantial power from a spinning BH, provided they are in the MAD state. The jet profile out to $$100\, GM/c^2$$ is roughly parabolic with a power-law index of k ≈ 0.43−0.53 during the sub-Eddington evolution. Both models show significant variability in the outgoing radiation which is likely associated with episodes of magnetic flux eruptions. The a* = 0.9 model shows semiregular variations with a period of $$\sim 2000\, GM/c^3$$ over the final $$\sim 10\, 000\, GM/c^3$$ of the simulation, which suggests that magnetic flux eruptions may be an important source of quasi-periodic variability. For the simulated accretion rates, the a* = 0 model is spinning up while the a* = 0.9 model is spinning down. Spinup–spindown equilibrium of the BH will likely be achieved at 0.5 < a*, eq < 0.6, assuming continuous accretion in the MAD state. 

Tidal disruption events (TDEs) around supermassive black holes (SMBHs) are a potential laboratory to study super-Eddington accretion disks and sometimes result in powerful jets or outflows which may shine in the radio and sub-millimeter bands. In this work, we modeled the thermal synchrotron emission of jets by general relativistic radiation magneto-hydrodynamics (GRRMHD) simulations of a BH accretion disk/jet system which assumed the TDE resulted in a magnetized accretion disk around a BH accreting at ∼12–25 times the Eddington accretion rate. Through synthetic observations with the Next Generation Event Horizon Telescope (ngEHT) and an image reconstruction analysis, we demonstrate that TDE jets may provide compelling targets within the context of the models explored in this work. In particular, we found that jets launched by a SANE super-Eddington disk around a spin a*=0.9 reach the ngEHT detection threshold at large distances (up to 100 Mpc in this work). A two-temperature plasma in the jet or weaker jets, such as a spin a*=0 model, requires a much closer distance, as we demonstrate detection at 10 Mpc for limiting cases of a*=0,R=1 or a*=0.9,R=20. We also demonstrate that TDE jets may appear as superluminal sources if the BH is rapidly rotating and the jet is viewed nearly face on. 

ABSTRACT We present general relativistic radiation magnetohydrodynamics (GRRMHD) simulations of super-Eddington accretion flows around supermassive black holes (SMBHs), which may apply to tidal disruption events (TDEs). We perform long duration ($$t\ge 81,200\, GM/c^3$$) simulations that achieve mass accretion rates ≳11 times the Eddington rate and produce thermal synchrotron spectra and images of their jets. Gas flowing beyond the funnel wall expands conically and drives a strong shock at the jet head while variable mass ejection and recollimation, along the jet axis, results in internal shocks and dissipation. Assuming the ion temperature (Ti) and electron temperature (Te) in the plasma are identical, the radio/submillimetre spectra peak at >100 GHz and the luminosity increases with BH spin, exceeding $$\sim 10^{41} \, \rm {erg\, s^{-1}}$$ in the brightest models. The emission is extremely sensitive to Ti/Te as some models show an order-of-magnitude decrease in the peak frequency and up to four orders-of-magnitude decline in their radio/submillimetre luminosity as Ti/Te approaches 20. Assuming a maximum VLBI baseline distance of 10 Gλ, 230 GHz images of Ti/Te = 1 models shows that the jet head may be bright enough for its motion to be captured with the EHT (ngEHT) at D ≲ 110 (180) Mpc at the 5σ significance level. Resolving emission from internal shocks requires D ≲ 45 Mpc for both the EHT or ngEHT. 

ABSTRACT We use the general relativistic radiation magnetohydrodynamics code KORAL to simulate the accretion disc formation resulting from the tidal disruption of a solar mass star around a supermassive black hole (BH) of mass 106 M⊙. We simulate the disruption of artificially more bound stars with orbital eccentricity e ≤ 0.99 (compared to the more realistic case of parabolic orbits with e = 1) on close orbits with impact parameter β ≥ 3. We use a novel method of injecting the tidal stream into the domain, and we begin the stream injection at the peak fallback rate in this study. For two simulations, we choose e = 0.99 and inject mass at a rate that is similar to parabolic TDEs. We find that the disc only becomes mildly circularized with eccentricity e ≈ 0.6 within the 3.5 d that we simulate. The rate of circularization is faster for pericenter radii that come closer to the BH. The emitted radiation is mildly super-Eddington with $$L_{\rm {bol}}\approx 3{-}5\, L_{\rm {Edd}}$$ and the photosphere is highly asymmetric with the photosphere being significantly closer to the inner accretion disc for viewing angles near pericenter. We find that soft X-ray radiation with Trad ≈ 3–5 × 105 K may be visible for chance viewing angles. Our simulations suggest that TDEs should be radiatively inefficient with η ≈ 0.009–0.014. 

ABSTRACT We present the results of nine simulations of radiatively inefficient magnetically arrested discs (MADs) across different values of the black hole spin parameter a*: −0.9, −0.7, −0.5, −0.3, 0, 0.3, 0.5, 0.7, and 0.9. Each simulation was run up to $$t \gtrsim 100\, 000\, GM/c^3$$ to ensure disc inflow equilibrium out to large radii. We find that the saturated magnetic flux level, and consequently also jet power, of MAD discs depends strongly on the black hole spin, confirming previous results. Prograde discs saturate at a much higher relative magnetic flux and have more powerful jets than their retrograde counterparts. MADs with spinning black holes naturally launch jets with generalized parabolic profiles whose widths vary as a power of distance from the black hole. For distances up to 100GM/c2, the power-law index is k ≈ 0.27–0.42. There is a strong correlation between the disc–jet geometry and the dimensionless magnetic flux, resulting in prograde systems displaying thinner equatorial accretion flows near the black hole and wider jets, compared to retrograde systems. Prograde and retrograde MADs also exhibit different trends in disc variability: accretion rate variability increases with increasing spin for a* > 0 and remains almost constant for a* ≲ 0, while magnetic flux variability shows the opposite trend. Jets in the MAD state remove more angular momentum from black holes than is accreted, effectively spinning down the black hole. If powerful jets from MAD systems in Nature are persistent, this loss of angular momentum will notably reduce the black hole spin over cosmic time. 

We explore the plasma matter content in the innermost accretion disk/jet in M87* as relevant for an enthusiastic search for the signatures of anti-matter in the next generation of the Event Horizon Telescope (ngEHT). We model the impact of non-zero positron-to-electron ratio using different emission models, including a constant electron to magnetic pressure (constant βe model) with a population of non-thermal electrons as well as an R-beta model populated with thermal electrons. In the former case, we pick a semi-analytic fit to the force-free region of a general relativistic magnetohydrodynamic (GRMHD) simulation, while in the latter case, we analyze the GRMHD simulations directly. In both cases, positrons are being added at the post-processing level. We generate polarized images and spectra for some of these models and find out that at the radio frequencies, both of the linear and the circular polarizations are enhanced with every pair added. On the contrary, we show that, at higher frequencies, a substantial positron fraction washes out the circular polarization. We report strong degeneracies between different emission models and the positron fraction, though our non-thermal models show more sensitivities to the pair fraction than the thermal models. We conclude that a large theoretical image library is indeed required to fully understand the trends probed in this study, and to place them in the context of a large set of parameters which also affect polarimetric images, such as magnetic field strength, black hole spin, and detailed aspects of the electron temperature and the distribution function. 

Aims.We investigated the polarization and Faraday properties of Messier 87 (M87) and seven other radio-loud active galactic nuclei (AGNs) atλ0.87 mm (345 GHz) using the Atacama Large Millimeter/submillimeter Array (ALMA). Our goal was to characterize the linear polarization (LP) fractions, measure Faraday rotation measures (RMs), and examine the magnetic field structures in the emission regions of these AGNs. Methods.We conducted full-polarization observations as part of the ALMA Band 7 very long baseline interferometry (VLBI) commissioning during the April 2021 Event Horizon Telescope (EHT) campaign. We analyzed the LP fractions and RMs to assess the nature of Faraday screens and magnetic fields in the submillimeter emission regions. Results.We find LP fractions between 1% and 17% and RMs exceeding 10⁵ rad m⁻², which are 1–2 orders of magnitude higher than typically observed at longer wavelengths (λ>3 mm). This suggests denser Faraday screens or stronger magnetic fields. Additionally, we present the first submillimeter polarized images of the M87 jet and the observed AGNs, revealing RM gradients and sign reversals in the M87 jet indicative of a kiloparsec-scale helical magnetic field structure. Conclusions.Our results provide essential constraints for calibrating, analyzing, and interpreting VLBI data from the EHT at 345 GHz, representing a critical step toward submillimeter VLBI imaging. 

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87^*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value ofτ = 0.08_−0.02^+0.03with a position angle ofξ = 50.1_−7.6^+6.2 degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87^*emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.